AV&1  9 //-** 


APPROXIMATE  FORMULAS  FOR  MANY-CENTER  INTEORALS 
IN  TilE  THEORY  OF  MOLECULES  AND  CRYSTALS* 

. 

uuuiyiui  j Molecular  Structure  uuU  Spectra 

Department  of  Physics 
The  University  of  Chicago 
Chicago  37,  Illinois 

PROBLEM  in  the  theory  of  molecules  and  crystals  is  the  evaluation  of 
nter  integrals  of  the  type 

( ab | cd ) - //  ^a(l)^b(l)[e2/r12]^c(2)V'd(2)dT1dT2  , (1) 

where  the  ^-functions  are  normalized  atomic  orbitals  associated  with  the  nuclei  a,  b, 

c,  and  d,  which  in  special  cases  also  may  coincide.  Useful  approximate  formulas  for 
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these  integrals  have  been  given  intuitively  by  Sklar  and  by  Mulliken  , who  have 
expressed  them  in  terms  of  overlap  integrals^ 

Aab  " / WT  " 6ab+Sab  <2> 

and  two-center  integrals  of  the  Coulomb  type.  Recently  Ruedenberg1*  has  shown  rigor- 
ously that  the  Mulliken  formula  may  be  considered  aB  being  the  first  term  in  an 
infinite  expansion  of  (1). 

At  the  Shelter  Island  Conference,  1951,  J.  E.  Mayer  pointed  out  that  charge 
densities  (expressed  in  atomic  units  with  e 1)  of  the  form  i>  * -A  .*  2 do  not  give 
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*This  work  was  assisted  by  the  Office  of  Naval  Research  under  Task  Order  IX  of 
Contract  N6ori-20  with  The  University  of  Chicago. 

^Permanent  address:  Institute  for  Mechanics  and  Mathematical  Physics,  The  University 

of  Uppsala,  Uppsala,  Sweden. 

1A.  L.  Sklar,  J.  Chem.  Phys . 7,  984  (1959);  A.  London,  J.  Chem.  Phys . 1£,  396  (1945). 

2R.  S.  Mulliken,  J.  chim.  phys . 46,  497,  675  (1949),  see  particularly  Eqs.  (63)  and 
(154b). 

^Ncte  the  difierent  meanings  of  the  symbol  S in  papers  by  different  authors;  we  have 

always  used  this  notation  for  the  devlat ion  from  normality  and  orthogonality  for  a 

set  of  functions,  and,  for  a **  b,  we  have  therefore  S » 0. 

aa 

**K.  Ruedenberg,  J.  Chem . Phys . 19,  1433  (1951). 
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any  contribution  to  the  total  oharge  and  that  these  "quasi-dipoles"  may  be  used  in  the 
interpretation  of  the  energy  expression^.  The  purpose  of  this  note  is  to  show  that  it 
is  possible  to  use  quantities  of  a similar  type  for  deriving  a generalization  of  the 
Sklar-Mulliken  formulas. 

For  the  sake  of  simplicity,  let  us  assume  that  the  oharge  distribution  has  a 
"oenter  of  gravity"  T,  whioh  is  situated  on  the  line  connecting  the  nuolei  a and  b 
with  respective  distances  a and  0 from  them,  and  that  similarly  ^c^d  has  a oenter  T' 
on  the  line  connecting  c and  d with  the  distances  y and  6 from  them.  We  will  then 
determine  the  coefficients  in  the  expression  X^^+X^2  from  the  condition  that  this 
quantity  should  have  the  same  total  charge  and  total  moment  as  the  density  whioh 

gives 

xa  " Aab(1+a/p)_1  ' xb  " Aba*1**/®)"1  (5) 

with  Aab  = Aba.  In  comparison  to  Mayer's  idea  this  means  that  we  are  here  considering 
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a quantity  y»a^b-A  ^ a _xb^b  whioh  doeB  not  contribute  to  the  total  oharge,  nor  to  the 
total  moment.  Neglecting  the  influence  of  higher  moments,  we  obtain  the  approximate 
formula 

Vb  " Va2+Vb2  » 

and  similarly 

Vd  " Vo2+Vd2  * (5) 

Substituting  these  expressions  into  (1),  we  get  finally 

(ab|od)  ■*  XaX0(aa|oo)+XaXd(aa|dd)+XbXc(bb|co)+Xbxd(bb|dd)  . (6) 

Hence  it  is  possible  to  evaluate  the  many-center  integral  (1)  approximately  in  terns 
of  the  overlap  integrals  Aab  and  Aod,  the  quotients  a/0  and  y/t>,  and  the  Coulomb 
integrals . 

For  symmetrical  oharge  distributions  with  a =>  0 and  y = 6 the  right-hand  side  of 
(6)  will  reduce  to  the  Mulliken  formula.  However,  in  the  case  of  unsymmetrioal  oharge 
distributions,  aB  for  instance  in  ionic  crystals,  formula  (6)  accr.?  to  be  an  essential 
improvement  of  the  previous  one,  and,  in  such  a case,  it  probably  also  indicates  a 

5See  also  0.  L.  Montet,  S.  P.  Keller,  and  J.  E.  Mayer,  J.  Chem.  Phys.  20,  1057 
(1952),  and  Keller's  and  Montet 's  Ph.  D.  Theses  (The  University  of  Chicago). 
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method  of  improving  the  convergence  of  the  Ruedenberg  prooedure  by  using  unsyrnmetrical 
multiplication  faotors. 

O 

In  order  to  derive  the  Sklar  formula,  we  observe  that  the  quantities  ^ab^T  and 

O 

Acd^T'  have  the  8ame  total  charge  and  the  same  total  moment  as  the  given  distributions 
^a^b  and  Vd*  respectively.  Negleoting  higher  moments,  we  get  therefore 

Vb  * Wt2  » Wd  " ^ed^T'2  1 <7> 

Here  V'.j  and  V'.pi  are  atomic  orbitals  associated  with  the  oonters  of  gravity  T and  T1, 
respectively,  whioh  have  to  be  chosen  in  a suitable  way1.  Substitution  of  (6)  into 
(1)  gives  then  the  approximate  formula 

(ab|cd)  - AabAQd(TT|T'T')  , (8) 

whioh  IS  essentially  the  Sklar  formula.  By  combining  the  forms  of  (4)  and  (7),  a 
still  better  result  may  sometimes  be  obtained  by  using  a three -parameter  formula  and 
a third  condition  related  toons  of  the  higher  moments. 

The  effect  of  the  higher  moments  in  general  and  the  correction  terms  in  (6)  and 
(8)  may  be  found  by  comparing  -these  formulas  with  the  exact  expressions  derived  in 
another  way®.  However,  a direct  comparison  with  the  numerical  values  of  the  integral 
(1)  in  special  oases  previously  obtained  by  us7  for  some  alkali  and  alkali  metals  and 

O 

with  Lundqv let's  material  for  LiH  shows  that  formulas  (6)  and  (8)  in  most  oases  have 
a surprisingly  high  aoouraoy,  put  also  that  there  exist  exceptional  Integrals  of 
certain  degenerate  types  whioh  are  better  treated  by  more  direct  methods. 

Formulas  (6)  and  (8)  are  intended  for  use,  e.g.,  in  the  MO-LCAO  method  based  on 
ordinary  atomic  orbitals  (AO).  If  the  molecular  orbitalB  instead  are  built  up  from 
orthonormal lzed  AO9,  the  use  of  (4)  and  (5)  for  the  ordinary  AO  will  lead  to  results 
which  are  related  to  those  reoently  intuitively  given  by  Parr10. 


®S.  0.  Lundqvist  and  P.-O.  Lowdln,  Arklv  FyBlk  ?.  147  (1951). 

7See  P.-O.  Lowdln,  J.  Chem.  Phya.  18,  365  (1950),  lg,  1579  (1951)  and  references  in 
these  papers , 

Q 

S.  0.  Lundqvist,  private  coiranunication. 

9P.-0.  Lowdln,  J.  Chem.  Phya ■ 18,  365  (1950). 

10R.  0.  Parr,  J.  Chem.  Phya.  20,  1499  (1952). 
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Finally,  we  wish  to  remark  that  the  author's  expressions11  for  the  total  and 
cohesive  energies  for  the  ground  state  of  a molecule  or  oryatal  having  its  total  wave 
funotion  approximated  by  a single  determinant  may  be  essentially  simplified  by  using 
formula  (6)  and  by  summing  the  contributions  from  all  matrix  elements  constituting  the 
coefficients  of  the  remaining  Coulomb  integrals;  in  this  process  it  is  even  possible 
to  treat  the  degenerate  Integrals  separately.  More  details  will  be  given  in  a survey 
of  the  various  methods  of  evaluating  many-center  integrals  to  be  published  elsewhere. 

I would  like  to  express  my  sinoere  gratitude  to  Professor  R.  S.  Mulliken  for  many 
valuable  discussions  and  for  his  great  hospitality  during  my  stay  in  Chicago. 


11P.-0.  Lowdln,  J.  Chem.  Phys.  1£,  1570  (1951),  Eqs.  (??)  and  (51). 
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